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ABSTRACT 

In this paper, we study the fermionic excitations near the quantum criticality using 
gauge/gravity duahty. This is motivated by exploring the non- Fermi hquid behavior under 
quantum phase transition. Fohowing the procedure of Liu, McGreevy and Vegh, we construct 
the gravity dual of "paramagnetic-nematic" phase transition in a continuum limit and study 
the non-Fermi liquid behavior under the quantum phase transition. We find that there exists 
non-Fermi liquid with the Fermi surface in the "paramagnetic" phase and the Fermi surface 
disappears in the "nematic" phase. 
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1 Introduction 



The traditional condensed matter physics is based on two corner-stones, namely Landau's 
Fermi liquid theory and Landau's symmetry-breaking theory. Landau Fermi liquids are con- 
trolled by a free Fermi gas fixed point with almost no relevant perturbations. While normal 
state of high-Tc superconducting cuprates and metals close to quantum criticality are exam- 
ples of non- Fermi liquids whose low energy properties are different from those predicted from 
landau's Fermi liquid theory (see [T]). In Landau's symmetry-breaking theory, thermal phase 
transitions are well developed based on so-called Landau-Ginsburg- Wilson paradigm. Quan- 
tum phase transitions are traditionally formulated within Landau paradigm which expresses 
the fluctuation of the order parameter in space and time directions. In recent developments 
[2] [3], new modes in phase transition which are inherently quantum and are beyond order- 
parameter fluctuations appear to emerge as part of quantum critical excitations. The proper 
description of the physics of strongly interacting fermions has long been a major issue for 
our understanding of condensed matter system near quantum criticality. The gauge/gravity 
duality has opened new avenues for studying strongly-coupled many body phenomena by re- 
lating a classical gravity in a weakly curved (d+l)-dimensional anti-de Sitter spacetime to a 
strongly-coupled d-dimensional quantum held theory living on its boundary. 

In this paper, we study the fermionic excitations near the quantum criticality using 
gauge/gravity duality. This is motivated by exploring the non- Fermi liquid behavior under 
quantum phase transition. Following the procedure of [1] [13!, ^® construct the gravity dual 
of "paramagnetic-nematic" phase transition in a continuum limit and study the non-Fermi 
liquid behavior under the quantum phase transition. At flnite temperature, we study the 
condensation of Yang-Mills fleld in the probe limit. At zero temperature limit, we consider 
the backreaction of Yang-Mills fields and find a fully gravity back-reacted solitonic solution 
to this system. Then we explore the non-Fermi liquid behavior in these backgrounds. We 
find that there exists non-Fermi liquid with the Fermi surface in the "paramagnetic" phase 
and the Fermi surface disappears in the "nematic" phase. 

An overview of the paper is as follows: In section 2, we set up the holographic model 
dual to paramagnetic-nematic phase transition. In section 3, we first study the system in 
probe limit to determine its phase structure and then we work out the gravity solutions 
which are dual to "paramagnetic" and "nematic" phases in zero temperature limit with the 
backreaction. One of them is the closed-form solution, and the other solution is numerical. 
In section 4, we analyze the Fermi surface behavior in different backgrounds. We conclude 
with a brief discussion of our results and future directions in section 5. In appendix A, we 
give some details about numerical calculation and show some diagrams. 

Note that when we use the word paramagnetic-nematic phase transition there is no sense 
in which the microscopic degrees of freedom of our our system consist of spins that are 
randomly oriented or aligned on a bipartite lattice. We only use the term to describe the 
symmetry breaking pattern. In the discussion of phase transition we are mainly interested 
in physical scenarios and dynamical mechanisms which arise from the holographic systems 
rather than detailed phase structure. The discussion of non-Fermi liquid behavior is similar, 
and we will not restrict to any specific theory. Since Einstein's gravity coupled to matter 
fields captures universal features of a large class of field theories with a gravity dual [15] , we 
work with this universal sector. Thus we take the so-called "bottom-up" approach, i.e. we 
just consider a certain type of field spectrum without referring to specific theory. We expect 
that if the gravity dual could be constructed top-down for such systems they would likely 
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contain ingredients similar to those in our description. 



2 Set-up of the holographic model 

In a recent paper [4j, the authors constructed a holographic model to study the quantum phase 
transition and they focused on realizing antiferromagnetic phase where the SU (2) symmetry 
is broken down to U{1) by the presence of a finite electric charge density. In many models of 
condensed matter physics, the spin rotations decouple from spatial rotations at long distances, 
and can be considered as global symmetry when discussing the low energy physics. In the 
holographic context, such a global symmetry is described by an SU{2) gauge field in bulk. In 
this paper, we study a holographic model to realize nematic phase in which the global SU (2) 
symmetry is broken as well as the rotation symmetry by the presence of current density. 
In the bulk, we consider SU{2) x U{1) gauge theory coupled to AdS gravity as the start 
point. The condensation of SU{2) gauge field breaks SU{2) gauge symmetry as well as the 
rotation symmetry. The bulk SU{2) gauge fields are dual to the background charge currents 
and charge density in the boundary. The U{1) gauge field is dual to the boundary chemical 
potential which does not play an essential role in our discussion. 

The bulk theory is described by gravity coupled to Yang-Mills field and Maxwell field with 
the Lagrangian 

= ^ + ^ - - \gmnG'^^^, (2.1) 

where TZ is Ricci scalar, R is the curvature radius of AdS4, F^^j^ is Yang- Mills field strength, 
and Gmn is Maxwell field strength satisfying 

Fmn = dMA% - dNA% + iq[AM, An]", 

Gmn = OmDn - QnDm- (2.2) 

We will restrict attention to the following ansatz for gauge fields. 

Am = A{r)T^dt + B{r)T^dx, (2.3) 
D,, = D{r)dt. (2.4) 

where q is the Yang-Mills coupling and r is defined as r"^ = ^ [a"" are the Pauli matrices), 
satisfying [t"", t^] = ze^^'^r^. It is straightforward to find the nonzero components of Yang- Mills 
field strength, 

FI = -FI = B'{t), F% = -FI=A{t), Fl = -Fl = qAB. (2.5) 

and the nonzero components of Maxwell field strength, 

Grt = -Gtr = D\r). (2.6) 

According to the anisotropy of Yang-Mills field ansatz in spatial direction, we choose the 
following ansatz for our metric, 

ds2 = -g{r)e-^^-'Ut^ + ^ + r'^ {c{rfdx'^ + dy^) . (2.7) 

gir) 
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With the above assumption, the equations of motion of Yang-Mills field can be reduced to 

\r I c J v^c^g 
B"+l<L-t_'L)B'+'^B = , (2.8) 

and the equation of motion of Maxwell field is given by 

D"+(- + ^ + -]d' = (2.9) 
\r 2 CI 



We use the definition of energy momentum tensor as 

1 '55'niatter 



T^. = - ^ 7 :: (2-10) 



and the nonzero components of energy-momentum tensor are given by 

Trr = ~A'' + j^B'' + -^{qABf-^D'^ (2.12) 
Ag Ar^c^ Ag^r^c^ 4g ' 

T,, = ]e^r^c^A'^ + ^B'^-^{qABf + ]e^r^c^D'^ (2.13) 
4 4 Ag 4 

T.. = \e-r^A--^,B- + £,iqABr + \e-r^D- (2.14) 

where T^x 7^ Tyy which is due to our anisotropic ansatz. Put the above result into the Einstein 
equations 

Rab - -^9abT^ - 3gab = Tab (2-15) 
and after a little algebraic manipulation the Einstein equations can be reduced to 

r c \g J r'^g'^c^ 

,/l c'\ fl 3c' c"\ ^,0 gB'^ e>^iqABf 

\r 2c J cr c J 4 4r^c^ Agr'^c'^ 4 



(2.16) 



It is to be noted that the above equations (2.16) have the following scaling invariance 

r — )• air, (t, x, y) — )• (t, x, y)/ai, g ^ afg, A ^ aiA, B ^ aiB, D ^ aiD. 
e^^a^e^, t 02*, A^A/a2, D ^ D/a2. 

x^x/a^, B^a^B. c ^ a-^c. (2-17) 

The second scaling symmetry can be used to set x = at infinity and the third one can be 
used to set c = 1 at infinity. Thus the metric solution has the asymptotic AdS4 form. 



4 



3 Solution to the holographic model 



3.1 Phase structure in the probe Hmit 



In this subsection, we will determine the phase structure of action (2.1) in the probe limit. 
The phase diagram of similar action with different boundary condition where one direction is 
compactified to a circle (and fermions are antiperiodic around this circle) has been discussed 
in a recent paper Compared their compactification model, the phase structure of our 

model is really simple. In order to study the phase structure analytically, we work in the 
probe limit where the Yang-Mills coupling is large and the backreaction of Yang-Mills on 
metric is suppressed. It retains most of interesting physics of the phase transition since the 



nonlinear interactions between Yang-Mills fields are retained. We start with Lagrangian (2.1 ), 
and then rescale Yang-Mills field as A]\j — )• ^^m- The Lagrangian (2.1) becomesj^ 



^ = ^ + J, - ^^PmnF^''''' - \GmnG'^^^ (3.1) 
where the definition of Yang-Mills field strength becomes 

i^Miv = dMA% - dNAlj + i[AM, A^r (3.2) 

In the large q limit, we can treat Yang-Mills as the probe and the background is determined 
by Einstein-Maxwell theory 

^em ~ ^ ~^ ~ '^^mnG^^ . (3.3) 

In the following, we will ignore the backreaction of Yang-Mills fields. It would be nice to 
include the backreaction, at least perturbatively to leading order. The Yang-Mills action is 

^ ' j4„ / — -77a TpMNa 



= --02 / d^x^gFtiMF'"''''''. (3.4) 



The analytic solution to Einstein-Maxwell theory is Reissner-Nordstrom-AdS black hole 

ds^ = -^dt^ + + -^{dx^ + dy^), (3.5) 



g{r) i?2 

with 



I)(r)=Ml-"^), (3.6) 



r 



^2 

g(r)=r^f{r), rj = — (3.7) 



f{r) = l-'\{l+v) + % (3.8) 



^In this subsection, we will restore the AdS radius R for dimension analysis. 
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The Yang-Mills fields in the Reissner-Nordstrom-AdS black hole background reduce to 

yl" + ^ A = 



The above Yang-Mills equations have the following scale invariance 

r^ar, {t, x,y) ^ {t, x,y) /ai, g ^ g, A ^ aA, B^aB. 



(3.9) 
(3.10) 

(3.11) 



Following the procedure of f8] we will study the condensation of B field analytically. By 
changing to the variable z = equations (3.9) and (3.10) become 



A" - B^A = 



with 



B" + ^2 Jr,z-s(i+v) B' + = 

l + 2*»7-2'^(l+?7) / '"r 



(3.12) 
(3.13) 

(3.14) 



and z = corresponds to the boundary and z = 1 corresponds to the horizon. Next we 
consider the boundary condition with these variables. Regularity at horizon requires 



A{1) = 0, B'{1) = 
and the asymptotic solution near the boundary reads 

A{z) = ao + aiz, B{z) = f3o + (3iz 



(3.15) 



(3.16) 



We fix the charge density rHCt2 = p and take f3i to be zero. Around horizon z = 1, we expand 
A and B as 



A{z) = A{1) + A'{l)iz - 1) + ^A"{l){z - If + ... 
B{z) = B{1) + B'{l){z - 1) + Ib"{1){z - If + ... 



(3.17) 
(3.18) 



Prom the boundary condition ( |3.15[ ), ^4(1) = and B'{1) = 0, and we set ^'(1) < and 

ed as 

(3.19) 



B{1) > for positivity of A and B. From the equation of motion A"{1) can be expressed as 

R'^B^ R^B{lfA'{l) 



A" 



z=l 



H 



giving 



A{z) = A {l){z - I) + 



R'^BilfA'il) 
24(ry-3) 



^-1 



(3.20) 
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The 2nd derivative of B can be calculated similarly as 

R"l ,2 4r?z-3(l + ??) ; 

-D 2 = 1 = ~^ A T7^. 2 = 1 



-B" 



R^B{l)A'{l)^ 



B\ 



2=1 



H 



2 = 1 



4(r/-3)2 



Thus we have 



4r|^(r/-3)2 

Around boundary 2 = 0, we expand A and B as 

A{z) = ao + ai2; + a2z'^ . 
B{z) = piz + p2z'^... 



(3.21) 
(3.22) 

(3.23) 



(3.24) 
(3.25) 



where we have used /3o = 0. Using equation of motion, the 2nd derivatives of A and B are 
given by 



02 = A"\ 



z=0 







/32 = S"|2=o = 
Then solutions around boundary become 

A{z) = ao + OLiz 
B{z) = /3iz 



(3.26) 
(3.27) 



(3.28) 
(3.29) 



where rnoii is fixed. In order to connect the solutions smoothly, we require the following 
conditions: 



/3i 



R^ba? 
16r|j(r)-3)2 



4r|,(r,-3)2 



Qq + jQi = — 2^ + 



8r^('?-3) 



ai = a 



2r^(^-3) 



(3.30) 
(3.31) 
(3.32) 
(3.33) 



where ^(1) = a and -B'(l) = b with (a < 0, 6 > 0). From (3.30) and (3.31) we can deduce 
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A = 3" 



From (3.32) and (3.33) we have 



2 ^ 2r|^(?7-3)(a-ai) 



(3.34) 



(3.35) 



From the AdS/CFT dictionary, the expectation value of dimension 2 operator is given by 

< O2 >= Piru/R^ (3.36) 
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Using the relation p = rnl^i and hawking temperature T = 4^^^ (3 — ij), < O2 > can be 
expressed as 

<0.>='^^#^4,/l-^ (3.37) 

where Tc is defined as 




V3p(3-r?) 

(3.38) 



8ttR 

We see that < O2 > is zero at T = Tc which is a critical point, and the condensation occurs at 
T < Tc- The critical temperature is the function of rj, so we can lower the critical temperature 
by adjusting the U{1) charge density. Tc — )• as ?? — )• 3, leading to a quantum critical point. 
Note that compared with [11], the phase structure of our model is simple. We have only two 
kinds of solutions the Reissner-Nordstrom-AdS black hole and a haired black solution. There 
is a phase transition between them and the critical temperature can be adjusted down to 
zero. 

Up to boundary counter terms, the free energy of the Yang-Mills theory is determined by 



the value of action (3.4) evaluated on-shell, which is the difference of free energy between the 



condensed phase and the uncondensed phsed. 

Ail — ^^condensed ^uncondensed — TS (3.39) 

where S is given by 

S = S'ymos + Sbdy + 'S'ct- (3.40) 

Using the equation of motion and the explicit form of metric, the on-shell action reduces to 

1 1 A^B^' 



/I 1 /" 

d^x --A'A\^=o + - J dz- 



f 



(3.41) 



The most efficient way to regulating S'ymos is working in the grand canonical ensemble where 
the chemical potential fi rather than the charge density is fixed. When working in the grand 
canonical ensemble, we do not need additional boundary terms and counter terms. Following 
the procedure of [9,j, we can calculate the AO using the above analytic approximation. 

1 2 1 . A^B^ 



An/V = --p' --J^ dz— (3.42) 

The analytic expression of AQ shows that the free energy of hair black hole solution is always 
lower than that of Reissner-Nordstrom-AdS black hole solution. The haired black hole is 
thermodynamically preferred. 



3.2 Zero temperature limit 

In this subsection, we will discuss the phase transition at zero temperature limit with back- 
reaction included. First we study the linearized perturbations of the extremal dyonic-like 
AdS black hole in order to analyze the instabilities of the backgrounds following the strategy 
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of [To] |12j . The general dyonic-like AdS black hole solution satisfying the constraint that 
g^yA^^A'^ is finite at horizon is given by 

c=l, X = B = ^ 

D = P. {i-D 

A = p2{l-l) 

(3.43) 



The horizon radius can be set to one using the first scaling symmetry of (2.17). The solution 



(3.43) is the most general background that preserves rotation symmetry. The temperature of 



dual field theory is given by 



47r 
Ibvr 

In the zero temperature limit, we have p'i + = 12. On the other hand, in order to study the 
gravity dual of "paramagnetic" phase, we must set /O2 = otherwise SU{2) symmetry will 
be broken. Thus in the background due to ordered phase, the Yang-Mills fields are inactive. 
The near horizon limit of the extremal solution reduces to AdS2 x M? with a metric, 

ds'' = _6(r - Ifdt^ + -r^—-^ + dx^ + dy^ (3.45) 

b[r — ly 

In the near horizon region, the equation of motion of B field reduces to 

djB + -dfB = (3.46) 
r 

where the effective mass of B field is ma^^ = and f = r — 1. In our case the Breitenlohner- 
Freedman bound for AdS2 is satisfied: 

m2g = 0>m|F = -^ (3.47) 

so there is no classical instability in this background. But this background bears the quantum 
instability due to the Schwinger pair production [5] which will lead to neutralisation of the 
black hole. In the dual field theory, this corresponds to the phase transition from ordered 
phase to the disordered one. 

When studied holographic superconductor, people have realized that the zero temperature 
black holes in Einstein-Maxwell-Higgs theory [7] [10] and Einstein- Yang-Mills theory [12] have 
the zero horizon size. Motivated by these facts, we assume the same is true for our model. 
Similar to that of |T0] [12], we choose the following ansatz near r — )• O"*" 

A^A^ir), Do{r) + Di{r), B^Bo + Bi{r), x~Xo + Xi(0, 

g^r'^ + gi{r), c~Co + ci(r). 

(3.48) 
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Terms with subscript 1 are subleading terms which go to zero faster than leading part. Sub- 
stitute the above ansatz (3.48) into equations of motion of Yang-MiUs fields (2.8), we have 



^2 / 2 Ai 



^^^.4o(r) ^ Ao(r) ~ e(-^\ a = qBo/co, 



(3.49) 



where we have used the constraint that Qf^uA^A'^ is finite at horizon. Following the similar 
procedures, to first order the equations (2.8) (2.9) and (2.16) have the following solutions 



8r2 



J^p-2a/r „ ^ ^2 _ eXOAgg 2a/r n ^ n 



(3.50) 



where we can set = Ij Xo = and cq = 1 by the scaling symmetries (2.17). After obtain the 



full solution, we can use the scaling symmetry again to set the asymptotical geometry to be 



AdS4. As a result, B component should vanish near the boundary. The gravity solution (3.50 ) 
not only breaks the C^(l)3 gauge symmetry generated by but also picks out x direction as 



special. So the background (3.50) captures the physics of the nematic phase which breaks the 



gauge symmetry as well as the rotation symmetry. 

Before concentrating on the fermionic probe, we can make some clarifications about gravity 
solutions here. If we want to find rotation invariance solution of our equations of motion 
K^i and KTO^, we must set c = 1 and 5 = for our ansatz p^l) and KTh. With 



these conditions, the equations of motion can be carried out with the closed-form solution 



(3.43). There are two free parameters pi and p2 in the solution. In this paper, we work 
in the zero temperature limit where pf + P2 = 12. First, when p2 = and p2 = 2\/3, the 
solution preserve SU{2) gauge symmetry and rotation symmetry around x and y directions. 
This background captures the physics of the paramagnetic phase. This background bears the 
quantum instability due to Schwinger pair production through there is no classical instability. 
Second, when /O2 7^ and pf + P2 = 12, this solution breaks the SU{2) symmetry down to 
U{1) but preserves the rotation symmetry. This background is dual to the ferromagnetic 
phase with external magnetic field nonvanshing since the external magnetic field is dual to the 
chemical potential of SU (2) gauge field. In this background, there are both classical instability 
and quantum instability. The classical instability comes from the violation of BF bound by 
effective mass of Yang-Mills field component B. This is consistent with the thermodynamics 
analysis [13| where the additional external field raise the free energy. Comparing with the 
above two solutions, our numerical solution is preferred at quantum level. In the following, 
we will focus on Fermi surface in the paramagnetic-nematic phase transition rather than in 
the paramagnetic-ferromagnetic phase transition. While the discussion about Fermi surface 
in the ferromagnetic-nematic phase transition is similar. 



4 Fermi surface behavior in different backgrounds 

Now we focus on the effect of the Yang-Mills field condensation on the non-Fermi liquid 
behavior. The basic idea, introduced in [11], is to introduce fermions in the above geometries. 
The Green's functions of fermions probe the existence and the properties of the fermi surface 
in the boundary theory. In our models, we suppose the fermions are charged under the U{\) 
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symmetry. The Yang-Mills fields affect the Fermi surface behavior only through background 
metric. 

To compute the spectral functions for fermionic operator in dual field theory, we need 
only quadratic action of ^jJ in geometry (3.43) and (3.48) 



spmor 



(4.1) 



where 



(4.2) 



and u^abM is the spin connection. We will use M, A^... and a,b... to denote bulk spacetime and 
tangent space indices respectively, and /x, z/... to denote indices along the boundary directions, 
i.e. M = {r,n). Although our discussion does not depend on dimension d, we will specify 
d = 3 which is dual to 2+1 dimensional field theory in the boundary. To analyze the Dirac 
equations following from (4.1), it is convenient to choose the following basis 

r'- = 



1 

-1 



jf" 
7^ 



(4.3) 



where ip± are the two-component spinors and 7^ are (2+1) dimensional gamma matrices. We 
can seperate the radial and boundary coordinate dependencies in ip as follows 



y± 

z± 



(4.4) 



First, we consider Dirac equation in the dyonic-like AdS black hole background (3.43) in 
the zero temperature limit since we are interested in quantum phase transition. The Dirac 
equation from the quadratic action is 



The background metric can be written as 

ds^-- 



(4.5) 



g{r)dt^ + lL_^+r\dx^+dy^), 



9{r) 



The Maxwell field D is given by 



1 + 



4r^ 



D = Dtdt, A = Pi 1 



Am = A{r)T^dt, A = p2{l 



(4.6) 



(4.7) 



The Yang-Mills field is 



In the zero temperature limit, pi and p2 satisfy the following condition: 

pI + pI = 12 



(4. 



(4.9) 
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Using the coordinate seperation (4.4), the Dirac equation becomes 




{dr =F rriy/g^) (t)± = ^iK^-^^cf)^, 



where 



K^i{r) = {-u{r), ki), u= . -^{u + - -)) . 

V ~9tt r 



(4.10) 



(4.11) 



We choose the fohowing basis 7*^ = ia2, 7"*^ = cii, 7^ = cxs for the 7^ and set A;2 = using 
the rotational symmetry of the dyonic-hke AdS black hole background. The Dirac equations 
can be reduced to the following two sets of decoupled equations 

-{dr =F = =Fi(fci - m)z=p, 

-{dr±m^/g^.)z^ = ±i{ki+u)y±, (4.12) 
We introduce the ratios ^+ = iy-/ Zj^, ^_ = — in terms of which the above equation 




(4.12) can be written as 




dri± = -2my/cHii± =F (fci Tu)± {ki ± u)^|. 



The retarded Green's function Gr is given in terms of the quantities by: 



Gn 

. 1 \ G22 



(4.13) 



(4.14) 



Prom now on we will drop the subscript 1 on momentum ki. The spinors ^± satisfy the 
in-falling boundary condition at the horizon which implies 



--1 = 1 ■ 



At w = 0, the in-falling boundary condition (4.15) should be replaced by 



m 



r=l,u=0 



k"^ + m? - 



Mi 



le 



(4.15) 



(4.16) 



With the boundary conditions, we can integrate equation (4.13) numerically to r — )• 00 



obtained the boundary correlation function. The imaginary part of the retarded function 



Gr is the function of k and uj for fixed m, qj and pi. The equation (4.13) which gives the 
location and property of the Fermi surface is nearly the same as that of Liu, McGreevy and 
Vegh |14j . So we find the similar Fermi surface behavior as that of Liu et.al. Therefore we will 
not repeat the numerical calculation here(we leave the numerical details in Appendix A), but 
the conclusion is obvious due to Liu et al.|14j. For some fixed w < 0, ImG22 has a sharp peak 
at some momentum value kp which implies the Fermi surface exist. Also from the scaling 
behavior of ImG22, it is different from that of landau Fermi liquid. 

Next we will consider the Dirac equation in the background whose near horizon limit gives 



(3.50). The key difference of this numerical background is that there is no rotation symmetry 
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in the x — y plane. Using the similar technique, the Dirac equation can be reduced to the 
following 



-Vgdr + m) 
{y/gdr + m) 



-—{-iLvy + + 

rc r 

{-luY) H \ 



V9 



rc 



Near r — )• oo, (j)± have the following asymptotic behavior, 



Br 



Cr 



,m~l 



+ Dr~ 



The coefficients D and A are related by a matrix 5, 

D = SA 

and then the retarded Green function Gr is given by [16J 

Gr = -i57° 



In the near horizon region, the equation (4.17) (4.18) becomes coupled equations for 



r{rdr — m)(l)^ 



-Z7 • 
-i7 • 



(4.17) 
(4.18) 

(4.19) 

(4.20) 

(4.21) 

(4.22) 
(4.23) 



where j ■ k = —lo^^ + kij^ + A:27^. These coupled equations can be solved using Bessel 
functions, and the solution satisfying the in-falling boundary condition at horizon is 



r ■^K,,. 



m+ 



a+ k'^ >0 
a_(_ oj > \k\ 
a_(_ uj < —\k\ 



(4.24) 



where a+ is an arbitary constant spinor. In the foUowings, we will focus on the spin averaged 
spectral function which is defined as 



1 



A{uj,k) = -Im(TrGR 



vr 



(4.25) 



From the numerical solution of equation (4.17) and (4.18) satisfying the proper boundary 
conditions, we can find the behavior of spectral function A{uj, k) with parameters m and q 
freely varied. For definiteness, let us now focus on m = and q = 1. In this case the spectral 
function A{uj,k) is a function w, ki and k2. We can plot the three-dimensional diagram of 
spectral function A{ijj, k) as the function of two variables of cj, k\ and ^2- Using mathematica 
we find that in the region — l<a;<l, — 3<A:i<3, and — 3 < A:2 < 3 the spectral function 
is regular which implies that the Fermi surface disappears. 

The disappearance of Fermi surface can be understood as follows. In the ordered phase, 
the Reissner-Nordstrom-AdS black hole sources an electric fiux towards infinity and there are 



13 



no charge- carrying fields in the spacetime. The fermion spectral function exhibit Fermi surface 
singularities. The origin of the Fermi surface is from the strong interaction between fermionic 
excitations and quantum critical bosonic modes in the IR. After the phase transition, the 
black hole loses their charge due to Schwinger pair production and the IR geometry AdS2 x 
M? collapses to the zero horizon gravity solution. In this paper, we do gravity calculation 
classically with including the backreaction of the pair production particles and so the charge 
density vanishes in our numerical background solution. If the backreaction is included, our 
numerical background may be modified to something like the electron star. It is interesting 
to study the fully back-reacted solution of disordered phase in our model. 

5 Summary and Discussion 

In this paper, we study the fermionic excitations near the quantum criticality in the holo- 
graphic model. We construct the gravity dual of "paramagnetic-nematic" phase transition in 
a continuum limit and study the non- Fermi liquid behavior under the quantum phase transi- 
tion. We also study the Yang-Mills condensation at finite temperature in the probe limit. We 
consider the Yang-Mills field in the background of extremal dyonic-like AdS4 black hole and 
find a fully gravity back-reacted solitonic solution to this system at zero temperature limit. 
The non-Fermi liquid behavior is studied in this background and it is found that the Fermi 
surface disappears under the quantum phase transition from gravity side. Though the Fermi 
surface disappears in the "nematic" phase, it is interesting to explore the non-Fermi liquid 
behavior further. The above discussion on Fermi surface behavior under quantum phase tran- 
sition can be generalized to other background such as Lifshitz case |17j . As mentioned earlier, 
it is interesting to construct our holographic model from the top down. In that case, we will 
know more about the field content and interaction form. 

In this paper, we construct the gravity dual of the nematic phase though the nematic 
phase. Study the Fermi surface behavior in some other kinds of quantum phase transition, 
such as the paramagnetic-ferromagnetic phase transition, is an interesting open question. 
In our model, the Yang-Mills fields obtain expectation value spontaneously which is not a 
standard problem in high energy physics since such an expectation value breaks the Lorentz 
invariance. One thing we should mentioned is that the gravity background dual to the nematic 
phase constructed in this paper is only a numerical solution and the analytical solutions for 
these kinds of backgrounds will be interesting. It would be desirable to study our construction 
at finite temperature with backreaction included. In that case the quantum phase transition 
becomes a thermal one and we can compare the holographic model result with the condensed 
matter result [18] . 
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Appendix A: Some Numerical Results 



In this appendix, we will explore the location of the Fermi surface in two different backgrounds 



(3.43) and (3.50) using numerical methods. First we solve equation (4.13) numerically with 



boundary condition (4.15) and (4.16) using mathematica NDsolve, and then we plot ImGa 
using formula (4.14) which gives the information about spectral function. For definiteness, 
we focus on m = 0, g = 1, and pi = and then our equation (4.13) is the same as that 
of Liu, McGreevy and Vegh [T3j. So we find the same diagram as that of Liu et al. (Figure 
1). The sharp peak in LTLG22 shows the singular behavior of spectral function which implies 





D.0 0.2 «.4 0.6 O.S l.D 



Figure 1: Left: Plots of ImGn as function of /c at cj = 0. Right: Plots of ImG22 as function 
of A; at u; = —0.001. A sharp peak in ImG22 is visible near kp ~ 0.92. 



the existence of Fermi surface. Secondly, about the numerical background which is dual to 




Figure 2: Left: Plots of ImGi2 as function of ki for fixed A;2 = and uj = —0.001. Mid- 
dle:Plots of ImG2i as function of ki for fixed A;2 = and w = —0.001 Right: Plots of ImGa 
(i=l,2) as function of ki and ^2 with uj scanning the region (-1,+1). 



"nematic" phase we need to NDsolve equation (4.17) and (4.18). Then using the definition of 
retarded Green's function, we can plot ImGu in the numerical background. We plot ImGa 
as function of ki and ^2 while oj scans some region. We find that in the region — 1 < u; < 1, 
—3 < fci < 3, and — 3 < A;2 < 3 the spectral function is regular (Figure 2) which implies that 
the Fermi surface disappears. While the off-diagonal components of spectral function have 
sharp peak which is consistent with the rotation symmetry breaking background. 
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